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A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 
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by 

PRAKASHKUMAR BABULAL POP AT 
Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 
April, 1986 

The present work is an investigation of the sheet metal 
deformation during blanking operation using finite element 
analysis. Blanking, punching, trimming, shaving and other 
operations for cutting metal sheets to shape using a die and a 
punch, forms an important group of press operations. These are 
the basic operations to several secondary operations such as 
drawing, coining, etc. 

Although blanking is one of the simplest and also one of 
the well studied press operations, there are a number of un- 
solved problems which are required to be taken up. While 
reviewing the available literature on this process, it is found 
that a large number of experimental and semi-experimental tech- 
niques are used in the study of the process. But the reported 
theoretical work is not exhaustive. The process of blanking is 
unsteady and a closed form analytical solution to this operation 
is not feasible. 

The finite element method is used in the present analysis 
because it is a versatile method and it is possible to predict 
the distribution of stress, strain, strain-rate etc. with suffi- 
cient accuracy using this technique. Based on the displacement 
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formulation, a general purpose elasto -plastic finite element 
programme is developed. The 'initial stress method' is used 
for the nonlinear analysis, with the material obeying von- 
Mises yield criterion. The parabolic isoparametric elements 
are used to discretize the continuum. 

The deformation during the process is assumed to be a two- 
dimensional plane strain one. Displacement-controlled analysis 
is performed by suitably prescribing the displacement increments. 
The material is assumed to be linear elastic-linear strain- 
hardening type. A large number of computer runs are taken by 
varying the process parameters. 

Punch-die clearance c , sheet thickness t, cutting conditions 

and the material properties such as elastic modulus E, Poisson's 

ratio Y» yield stress strain hardening parameter H' and 

yp 

local fracture strain are considered as the important para- 
meters of the process. In thi« work the effects of the fra- 
ctional punch-die clearance c/t, t, H' and on the punch 
force, punch penetration (at crack initiation), deformation 
pattern, optimum clearance etc. are investigated. The punch 
and the die surfaces are assumed to be flat, smooth and possess 
sharp edges. 

A crack is introduced into the sheet by using the follow- 
ing crack initiation criterion. This states that, a crack is 
initiated at any point in the material where the effective 
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strain becomes equal to or greater than the specified local 
fracture strain of the material. It is also postulated that, the 
crack propagation direction coincides with the direction of the 
maximum shear strain at tho crock initiating point. The execu- 
tion of the programme is terminated as soon as an indication of 
the crack initiation is obtained. The results, thus obtained, 
are used to plot a number of graphs showing the effects of the 
various parameters, the deformation patterns, the strain dis- 
tributions etc. From these an attempt has been made to find 
out tho optimum punch-dio clearance, maximum punch force, maxi- 
mum punch penetration etc. for different combinations of the 
process parameters. 


Based on this analysis it can be inferred that the maximum 
punch force is almost independent of the punch-die clearance. '■ 
With larger value of H’ , the punch force required (for the same 
punch penetration) increases, and finally to initiate a crack 
more punch force is needed. The variation of sheet thickness 
is also found to yield similar results- But it is observed that 
variation of t and H* has very little effect on the percentage 
punch penetration for a given value of c/t ratio. The percen- 
tage punch penetration (at crack initiation) increases with the 
increase in c/t ratio and the ductility (e^) of the material. 

It is found that the above observations agree well with the 
experimental results obtained by previous researchers in this 
field . 
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When the clearance between the punch and the die is optimum, 
the cracks originated from the punch and the die corners meet_ 
each other after propagating along a straight line* resulting in 
a clean edge. As it was not possible to get crack propagation 
path from the finite element analysis the following presumption 
is made to get the optimum clearance. It states that when the 
clearance is optimum, the line passing thi'ough the corners of 
the punch and the die (at crock initiation) coincides with the 
direction of the maximum sliear strain. By making use of this 
criterion optimum punch-die clearance for different values of 
ductility arc obtained. From these an estimate of the maximum 
punch penetration for each value of the optimum clearance is also 
made. The results are examined critically and verified with the 
available experimental and theoretical results. It is found that 
they are in close agreement with the experimental results and 
the trends match with the previously obtained theoretical 
results. The optimum clearance reduces as the ductility of the 
metal is increased. .So it can be concluded that the more ductile 
metals, e.g., Aluminium require less punch-die clearance com- 
pared to the less ductile metals like mild steel, for satisfa- 
ctory blanking. The optimum percentage punch penetration does 
not vary much with the ductility of the metal, and its value is 
found to remain around 30 /,, as obtained in previous experimental 


works 
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Finally, to visualise the deformation mode of the material 
in and around the clearance gap, the patterns of the nodal dis- 
placements at each incremental step are plotted. It is in- 
fered that the material near the punch and the die corners 
deforms most severely. This clearly brings out that the cracks 
should initiate near the punch and the die corners as observed 
in practice. Similar conclusioh is also drawn from the strain 
distribution contours which are plotted in the clearance gap, 
i.e., the maximum value of strains occur near the punch and the 
die corners for all cases considered. 

At the end it is concluded that with the help of this 
technique, it is possible to predict with sufficient accuracy 
the values of the optimum punch-die clearance and the press 
load required, without performing the experiments. Although 
attempts have been made by other researchers to analyze theore- 
tically the blanking process, but the results obtained by the 
finite element analysis, proposed in this work, are found to be 
much more accurate and the trends obtained regarding the effects 
of the various parameters using this technique are in closer 
agreement with the experimental work, when compared with pre- 
vious theoretical work. 



CHAPTER 1 


INTRODUCTION 


1.1 INTRODUCTION 

The technology of sheet-metal press-work emerged with the 
development of steel industry, and to a large degree we owe 
our present standard of living to the production of stamped 
metal parts. A large number of metal parts are manufactured 
on press, which are being used in household goods, automobiles 
and other items. 

Most of the operations for shaping sheet metal into 
finished products can be described as either cutting or forming 
operations. These are generally referred to as stamping, press 
working or punch press working of sheet metal. Such processes 
are very economical as the desired shape, size and finish can 
bo obtained without any significant loss of material. Moreover 
a part of the input energy is fruitfully utilised in improving 
the strength of the product through strain hardening. 

Blanking, punching, trimming, shaving and other operations 
of cutting metal sheets to shape using die and punch, forms a 
distinct group of press operations, the subdivision of which 
have much in common. This group is the largest among all the 
groups of press operations, the simplest and also the best 
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understood. The above mentioned operations are basic to 
several secondary operations such as drawing, coining etc. 

Still there are a number of unsolved problems and several 
complications which requires attention [1,2]. 

1.2 BASIC MECHANICS OF BLANKING 

In punching and blanking processes [3], a finite volume 
from a sheet metal is removed by using a die and a punch. 

The shape and the size of the portion removed are determined 
by the geometry of the die and the punch. If the final pro- 
duct happens to be the removed portion, then the operation is 
termed as blanking; on the other hand, if the pierced sheet 
metal is the final product then the operation is called punch- 
ing. However, the basic mechanics of metal removal is same in 
both the operations. 

Fig. 1.1 shows a simple punching operation. The job is 
hold by job holders to prevent any distortion and to provide 
a rigid support. The punch and the die corners are sharp, i.e., 
without any radius, as the objective in this process is to 
achieve rupture of the material. A clearance c is provided 
between the punch and the die. Hence the die diameter 
= dp+2c where d^ is the diameter of the punch. 

The process of sheet metal blanking is studied with two 
different aims. In one case the emphasis is purely on the 
study of the mechanics of the deformation process, while in the 
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Fig.I.I Details of punching process 
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other the effects of the various process parameters are 
studied. The most commonly accepted theory for blanking 
operation is based on the concept of shear failure of the 
material lying between the punch and the die. This theory 
results in a very simple formula to determine press load, 
which is given as follows : 

P=LxtxS (1 •la) 

max s 

or 

P ^=TTXdxtxS^ (1.1b) 

msix s 

where 

p^ax ~ maximum punch force to shear any shape with 
a sharp flat punch and a die having proper 
clearance between them. 

L = Length of the cut, t = thickness of the metal » 

d = Diameter of the blank, 

S = Ultimate shearing strength of the metal, 
s 

In this theory it is not possible to incorporate the 
effects of the other process parameters such as the amount of 
clearance between the die and the punch etc. 

Fig. 1.2 shows the nature of metal deformation as the 
punch penetrates into the workpiece. As the downward motion 
of the punch progresses, the workpiece material deforms and 
is pulled down by the punch movement causing the grains to 
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elongate near the punch corner B [S]. A similar type of 
deformation takes place near the other die corner D. When 
the grain elongation or the local natural strain in the 
surface fibre AB (of the workpiece) reaches a limiting value, 
the fibre ruptures. Since the local strain is maximum at 
the corner, a crack opens up at the punch corner* After 
this, with a slight movement of the punch, the inner fibres 
also got ruptured. Thus the fracture line BY propagates, 
following a path along which the successive inner fibres 
attain the fracture strain. A similar crack also propagates 
from the corner D (because of the symmetry of the deformation 
geometry) . Now if the amount of clearance c is optimum, then 
the two fracture lines meet and a clean edge is obtained after 
the operation (Fig. 1.3a). If the clearance is too small, 
then the fracture lines miss each other and the secondary 
deformation takes place, resulting in an unclean edge 
(Fig. 1.3b). Fig. 1.3c shows what happens when the amount 
of clearance is too large, a significant amount of drawing 
action takes place and the quality of the workpiece is again 
quite poor. 

1.3 PREVIOUS WORK 

The first scientific paper on blanking appears to be that 
by Anthony [4] published in 1911. It describes a pressure 
recording apparatus suitable for use with a punching machine. 
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Fig. 1.2 Deformation of workpiece with punch travel 



(a) (b) (c) 

(a) Correct clearance , (b) Too little clearance 
(c) Too large clearance 


Fig. 1.3 Fracture lines with different clearances 
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The principle of Emery testing machine was used in the design 
of the said apparatus. The apparatus was designed to enable 
an indicator diagram to be taken under the conditions of actual 
shop practice by applying the mechanism to punching machinery 
which were used during those days. The indicator being hydrau- 
lically operated, there was possibility of reduction of pre- 
ssure due to leakage if sufficient care was not taken. It is 
also a wellknown fact that the indicators, used under such 
conditions, are most likely to leak or stick and either con- 
dition would affect the results. Of course, it is mentioned 
that preventive action was taken by frequent calibration, which 
removed any reduction in pressure due to leakage. It is added 
that by the continued tapping of the indicator cylinder, irre- 
gularity due to sticking of piston was almost entirely elimi- 
nated and resulted in a high degree of precision. Such hydrau- 
lic fcuru- lever type indicators are only suitable for slow speed 
operation. Chang and Swift [s] have reported a work in which 
stresses, metal flow pattern, and crack propagation during the 
shearing of metal bars of most industrial metals is studied by 
performing experiments. The shearing of half inch thick metal 
bars of different ductility were performed. The bars were 
sheared by the process of double shear which resembled double 
cropping operation. Some of the variables which were varied 
during the operation are (i) clearance between the die and the 
punch, (ii) the radius of cutting edges ground on the die and 
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the punch (i.e., dullness of the cutting edges), (iii) the 
amount of prior tensile straining given to the bar before 
experiment etc. The shearing load and the punch penetration 
values, were obtained by means of a pressure gauge and a dial 
indicator, respectively. Fig, 1.4 shows the typical load vs. 
punch penetration curves obtained for various metals, from the 
experiments. It is mentioned further that the development of 
strain and progress of fracture were studied by subjecting 
several specimens of each metal, to a different punch pene- 
tration depths, i.e,, one specimen was sheared to the point 
of maximum load, while the others wore sheared to the stages 
intermediate between the point of maximum load and complete 
fracture. By careful measurement of the deformed grid pattern 
with the help of a travelling microscope, the distribution of 
strain was obtained in each case. At the end, following 
important conclusions wore made. There are two different modes 
of severence in a normal shearing operation, (i) ideally 
ductile materials suffer bodily sliding action along the 
line of shear and final severence occurs at full penetration 
and (ii) materials which arc not ideally plastic develop cracks 
of tensile origin which tend to propagate at a small angle to 
the line of shear. After this Chang [6] has done the experi- 
ments on the shearing of metal blanks. In this work the 
results of investigation on the effects of clearance, tool 
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Fig. 1.4 Load penetration curves 
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shear and tool curvature in the shearing of circular blanks 
is presented. The thickness of the various materials tested 
were in the range of 0.036 to 0.5 inches. The variations in 
the diameter of the tools employed were of the order of one 
to six inches and the clearance between the punch and the dies 
were varied between 0 to 58 percentage of the material 
thickness. The conclusion which was reported is that the 
basic modes of fracture in blanking with circular tools are 
essentially the same as those found in the bar shearing 
experiments. One of the observations made during the experi- 
ments was that as the diameter of the tool was increased, 
there was a systematic decrease in the nominal ultimate 
shearing stress. This indicates that in the range of the 
clearances used in industry, t he sheari ng stran.ath of the 
material blanked is an inverse_ function of the tool diameter. 
Another observation that was made is the nominal ultimate 
shearing strength of all the metals tested except lead de- 
creased with the increase of clearance. This effect was most 
pronounced with the clearance less than 10^. For lead the 
results obtained were rather uncertain. The work done in 
blanking of mild steel decreased with the increase of clearance 
initially and then remained steady or even increased with fur- 
ther increase in clearance. A clearance in the range of C^ito 
10^ is recommended in order to obtain a clean fracture with 
minimum work required to be done, while sheet blanking with 
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circular tools. It is further stated that the work required 
for shearing increases with the increase in the curvature of 
the tool. Also the fall in the maximum blanking load with 
increase. In the angle of shear was found to be almost linear 
and recorded a reduction of about 75^ in maximum load. Although 
the shear on the punch increases the punch penetration at maxi- 
mum load and at failure, but it helps in smoothing out the 
operation by eliminating shocks. Crane [2] has given a detailed 
discussion on the effects of the various parameters on the 
shearing process, which are similar to those observed by the 
previous workers mentioned above. It is mentioned that the 
shearing process in a punch and a die resembles testing the 
specimen in tension very closely. But it is stated that any 
definite relationship between stresses in blanking and tensile 
test was not found. Biegel [7,8] has studied the effects of 
variation in the clearance between the punch and the die on 
the forces and work required during stamping operation. Expe- 
riments were performed using a single circular punch (without 
providing any shear) in combination with seven different dies. 
By mounting the punch to a strain gauge type dynamometer the 
force— time trace was obtained on an oscilloscope and the same 
was recorded by a Polaroid Land Camera. This measuring method 
is quite superior to that using mechanical cum hydraulic auto- 
graphic recorders. Finally it was concluded that the maximum 
shearing force is a function of the area of the material being 
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sheared, the material strength and is independent of clearance. 
But the amount of work required for shearing is a function of 
the shear area, material strength and also the clearance. At 
the end it is mentioned that to reduce the wear and tear of the 
punch and the die, the tool material must have high strength, 
good wear resistance and should be completely free of imper- 
fections in the region of cutting edges. Hugo [9] describes 
the stamping die performance and quality levels of parts pro- 
duced. In general there are two methods of measuring stamping 
die performance viz., one based on the number of parts produced 
per press run and the other based on the number of parts pro- 
duced per sharpening grind. Five classes of cut edge characteri- 
stics were recommended for measuring the quality levels of the 
parts produced (Table 1.1). These five classes viz., A,B,C,D 
and E can be distinguished in terms of the type of the burr 
produced (i.e., tensile, compressive, extrusion or combined), 
and the amount of roll over (i*e*, minimum, moderate, large or 
very large). It is added as the punch-die clearance is in- 
creased the part quality level shifts from type A to type E. 

Also by limiting the productivity measurements in terms of 
burr height for each class the individual effects of the many 
contributing factors, in the stamping operation can be properly 
evaluated. Although this is a good technique to assess the 
quality of stamped metal parts, but it is difficult to measure 
burr height of each and every stamping in mass manufacture. 
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Table 1.1 


Five Cutting Edges for Measuring Part Quality Levels 

Punch-to-Die Clearances 


Cuttinfi Edftts 

JOe to Joe Comer 
irtiktiown 

Burriutd Etimlehed 

A 

Burr: Tensile 
and Compressive 

Min. 

ftoiV^ver 

r 

B 

Burr: Tensile 
and Cenvossive 

Min. 

ReMver 

-T — 1- ■ 

C 

Bmt: Tensile 

Moderite 

RotLOver 

0 

Burr: Tenalie 

Litm 

Retl-Over 

T ■ I-. 

E 

Burr: Tenatie 
and Extrusion 

Very Large 

RoMvw 

. ^ 

Imroisee as Cutting 

Eitges Break Down 

Dtanrtlnft ttf Run 

eee- 

\\ 

fewii'I'i'.V 



ininUm mMJfWl VI rtl#fl 

Linger Runs as "" 

Cieahmce increaeos 

bsx4c 

Doubie 

Shear 

Secondary 

Shear 

General 

Purpoee Hole 



Ejector Punches Required 

Materia! 


Total Cteerance in Portent of Materiai TMcknees 


Aluminum-^hard 

1-2 

KaH 

18-20 

25-26 

40 pet max. 

—BOft 

1- 2 

■oh 

12-16 

16-20 

34 pet max. 

Brass— % hard 

1-^ 2 


12-16 

18-22 

48 pet max. 

— 'annealed 

3 


12-16 

16-20 

42 pet max. 

Bronze, phosphor 

3~ 5 

7-10 

20-24 

25-27 

50 pet max. 

Copper— % hard 

2~ 4 

e-10 

12.-16 

16-22 

50 pet max. 

— tfinealed 

1 1- 2 

4- 8 

10-14 

16-18 

50 pet max. 

Lead 

S- 5 

8-12 

16-15 

16-20 

44 pet max. 

Magnesium 

1- 2 

3- 5 

7- 9 

10-14 

32 pet max. 

Stainless steel 

2- 4 

8-10 

18-22 

25-27 

46 pet max. 

Steel, CBS (1020) 

2 - 4 

10-14 

16-20 

23-25 

42 pet max. 

high carbon 


22-26 

26-32 

34-38 

! 50 pet max. 
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Some researchers have applied visioplasticity technique 
of analysis during blanking operation. This is a good semi- 
experimental method to visualise the deformation of the metal 
sheet as the punch penetrates into the sheet metal. The for- 
raability of metal can also be judged from such an analysis. 
Keeler [lO,ll] has used circular grid system to evaluate 
sheet metal formability. It is claimed that circular grids 
permit immediate and simple measurement of the maximum elonga- 
tion of the sheet at any location. The grids were marked on 
the sheet by electrochemical etching. As the material was 
strained, a circle became an ellipse whoso major and minor 
axes were measured. Knowing the diameter of the original 
circle the strains were calculated. Santos and Organ [l2] in 
the beginning of the analysis have given a brief review of 
the principles involved in determining strains and stresses 
during the observation of the deforming grid. The equations 
for relating incremental deformations of grid, in two dimen- 
sions of an originally orthogonal cartesian grid, to the incre- 
mental plane strain are stated and explained. Two approaches 
are described for calculating the strain increments from the 
measurement of distorted grid as follows : 

(i) The incremental displacement U, at a grid point is de- 
duced. This U is resolved into components u and v which are 
parallel to the x and y axis respectively, from which the cur- 
ves of variation of u and v w.r.to x and y both are plotted. 
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Then the resulting slopes a^re used to 

obtain infinitesimal strains in the x and y directions, i.e., 

*5 Sy and 6y^y» This method involves graphical or numeri- 

cal differentiation and therefore demands the application of 
the tedious process of smoothing . 

(ii) In the second method the local strain components are 
determined directly by noting the relative distortions along 
convenient axes. Then the direction cosines of these axes 
with respect to the x and y axes are determined. Finally the 
strains are transformed w.r. to the x and y. 

This second method was used for calculation of strains in 
the work reported. It is claimed that there was no previous 
study of the ductile fracture processes using the visioplasti- 
city method, although attempts wore made only to give the qua- 
litative idea of the distribution of deformation patterns. 

Fig. 1.5 and Fig. 1.6 show the proposed slip line fields for 
initial and lateral stages of the process respectively. It was 
found that deformation mechanism implied in the resulting slip- 
line field was in accordance with the published observations on 
the behaviour of the material, subjected to that process. 

Kasuga and others [13,14] reported application of visioplasti- 
city method combined with strain gauge technique during blank- 
ing of circular specimens. The grid was photoeched on one 
half of the specimen which was in two halves, and the strain 
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gauges were applied on the upper and the lower surfaces. 

This technique made it possible to get the flow pattern of 
the metal through out the whole specimen, and even in the 
elastic region. A technique of the virtual displacement of 
punch was also applied for the analysis of metal flow. From 
tho results of this analysis it is reported that as the 
clearance becomes larger, the danger of crack initiation at 
less penetration increases. A simplified shearing model with 
scissor like tools was also proposed. The general drawback of 
using the visioplasticity method is that the process is to be 
interrupted after small penetration of the punch. Also blanking 
being nonsteady state process in which metal is continuously 
strained, the visioplasticity techni^que may not reflect the 
correct behaviour of the metal during the process , 

The critical review of the progress of research in the area 
of slow and fast blanking at ambient and high temperatures is 
presented by Johnson and Slater [15] . Only the literature 
available in English was reviewed. It is claimed that there 
were relatively few papers which add significantly to tho under- 
standing of the effects of the parameters involved in the blank- 
ing process. It is mentioned that modern blanking practice aims 
at producing the finished products which do not require any fur- 
ther machining or finishing-. There is a demand for high rate of 
production together with a high degree of accuracy, good edge 
finish and consistency in production. Only experiments under 
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laboratory conditions are performed at extremely slow or quasi- 
static conditions which give rise to a mean shear strain rate of 
the order of lo"* to lo'^ s'^- FiS- =hows a schematic repre 

^ ^ ^ ^nnch displacement diagram. It is said 

sentation of punch f orce-p^Ji^^" 

X ) ^ 1 ^ -From the punch and the die corners 

that cracks are developed 

XX ui 1 • X- -r-cpiches a maximum value and they pro- 

after the blanking force reac-n-^. 

X X j u xu^^ After performing the dimensional 

pagate towards each other. ^ ^ 

analysis of the parameters involved in circular blanking at 
constant strain rate the following expression is given : 


p^ax/"’’ dt = 0[(c/t) 


( 1 . 2 ) 


where 


max 


d = diameter of the blank 
t = thickness of sheet metal 
c = punch-die radial clearance 
T = current shear stress of the material 
= maximum blanking f°tce required . 


XX - xxu xx^xu-f- r^n. (1‘2) may be used to study the 

It IS further stated that eq» v-*- / / 

, , X • XI -,r.irlnq force and energy required on 

dependence of maximum blanKii‘y 

various parameters. Following conclusions were made from a 
number of experiments performed and review of literature : 

/ . X , . . . x X.; r- blanking force is independent of the 

(i) Maximum quasi-static Dirf“ y 

T T . X +v,c:>xp is a rising characteristic as the 

radial clearance but there j-p « a 

. j- 1 1 approaches zero, (ii) With the 

percentage radial clearance 
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a> 
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Maximum punch force 


Initiation of crack propagation 


Secondary cracks probably formed 



Frictional resistance between 
(a) Blank and hole 
Cb) Punch and hole 
(c) Blank and die 
Frictional resistance mainly 
between blank and hole 

Separation of Wank from shee 


Frictional resistance mainly 
between punch and hole 


Punch displacemer 


ig.1.7 Schematic representation of punch force -punch 
displacement autographic diagram 
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increase in the ratio of tool diameter to material thickness, 
the maximum blanking force per unit sheared area decreases. 
Regarding the effect of speed on the process it is stated that 
generally greater energy is required for dynamic blanking in 
comparison to quasi-static process. This can be attributed 
mainly to the effect of strain-rate on the yield shear stress 
of the material, if it is assumed that the mode of plastic 
deformation and the fracture are independent of speed. As 
regards the phenomenological features of the blanking process 
it is commented that the blanking of closed contour is a com- 
plex non steady-state process involving a triaxial stress 
system. It is well known that the process of producing blank 
with a die and a punch involves both plastic deformation and 
fracture? out of the two the former is reasonably well under- 
stood, but not the latter. It was observed that by providing 
small radius to the tool, the effect on the quasi-static 
blanking force was almost negligible. Further it is mentioned 
that the non-dimensional edge taper is approximately directly 
proportional to the percentage radial clearance. When the 
blanking of mild steel is carried out at high speed it was 
observed that there was a considerable reduction in distortion 
of the blanked piece in terms of ’doming’ and ’dishing'. In 
the later half of the reported work, the results of the axi- 
syrametric blanking at elevated temperatures is described. 
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Tatsuo and others [16,17] have reported the work of circular 
punching, two sided cutting and blanking using profiled tools, 
which was undertaken mainly to investigate the effect of profile 
on the resistance to shear and the quality of the sheared sur- 
face obtained. It is concluded that the material under the 
punch is deformed by indentation of punch in the early stages 
of the process, then a sort of forward-backward extrusion and 
ultimately by the usual process in punching of small holes. 

Blanking, bending and deep drawing are some of the most 
important techniques in press working. Among these techniques 
the plastic deformation of the sheet material in the process of 
bending and deep drawing has been studied well. However, for 
blanking process, where the material in the narrow region bet- 
ween the edges of the die and the punch flow unsteadily, litera- 
ture on theoretical analysis available is not adequate. The 
theoretical work done on blanking by Maoda [21] is reported by 
Jimma [l8,19] and Masuda and Jimma [20] together with their 
own work. Jimma applied theory of limit analysis for the theore- 
tical investigation of the blanking process. The process of 
circular blanking was considered, and the theory of plane strain 
deformation was used, assuming the material to bo a non- 
hardening plastic rigid body. It is reported that the process 
was analysed in two parts. In part I of the analysis following 
work is described : 
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(i) The sheet material was separated into several regions of 
constant stress by using the straight lines of discontinuity, 
and the statically admissible field was analysed, (ii) The 
lower limit of blanking force against each punch penetration 
was calculated, (iii) An estimate of the development of crack 
in blanking process is explained. 

In part II of the analysis, the plastic deformation of the 
sheet material in the blanking process was divided into four 
stages, viz., the first stage (obtuse plane or shear drop), 
the second stage (bright plane), the third stage (depre- 

ssed plane) and the fourth stage where the separation due to 
slip off occurs in the material, after the complete penetra- 
tion of the punch into the sheet. All these four modes of 
deformation for ideally ductile metal is shown in Fig. 1.8. 
Thereafter the kinematically admissible velocity field in the 
material was analysed for each stage and the upper limit of 
blanking force against punch penetration was calculated. Follow- 
ing hypothetical mechanism is also mentioned regarding the ini- 
tiation and growth of the microcrack during the process. Micro- 
cracks occur first in the material near the tool edges where the 
strain hardenability of the material become too low after severe 
distortion. When the microcrack moves from the tool tip to the 

i, 

side of the tool the normal stress across the plane of the crack 
becomes tensile and, thus, microcrack opens up. The conditions 
put forward for a microscopic crack to develop are : (i) the 
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(a) The first stage 


(b) The second stage 


Obtuse plane (Shear droop) 



(c) The third stage 


(d) The fourth stage 



Fig. (.8 Four modes of deformation in the shearing 
processes of ideally ductile metals 
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shearing strain along the slip plane near tool edges should 
reach a limiting value at rupture, and (ii) the normal stress 
across the plane of raicrocrack should become zero. Finally 
a number of experiments were performed to verify the theoreti- 
cal predictions and results which were obtained are reported 
to be quite satisfactory. 

It is observed from the literature reviewed so far that no 
one has put forward any explicit mathematical relation connect- 
ing the material properties, sheet thickness, and clearance, 
with the maximum blanking force, etc. Ghosh and others [22,23] 
have put forward a simple model for the deformation of the metal 
in the clearance gap during the blanking operation, from which 
a mathematical expression connecting the variables of the pro- 
cess has been tried. It is mentioned that the metal fibres rup- 
ture when the local strain in the material reaches a critical 
value, and using this concept it is claimed that a theoretical 
relationship between the optimum clearance and the other para- 
meters has been developed. " . 

Some parts of the said analysis is presented here. The 
relation between the active strain, and the locally obser- 
ved strains and ^2 was asVumecl'~t''b be related through the 

following expression [24], 

1-^1 B 2 I) 


s'l = V3( e 


(1.3) 
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This was simplified for the case of blanking, and the 
criteria for fracture was taken as 


or 


= 


” ®f 


( 1 «3a) 

(1.3b) 


where e^ is the engineering fracture strain of the material 
blanked. Fig. 1.9 shows the details of the assumed deforma- 
tion pattern for optimum clearance, when fracture has just 
started at point B. From the said figure and assuming linear 
variation of strains along the fibres it was shown that 

oA 

eg = ^ -2 + (1.4) 

o 

where Gg is the observed engineering strain at B and c^ is the 
optimum clearance. 

Similarly the expression for strain at S was obtained as 
% -2 + (2A /cq) 

= (1.5) 

1 + (2A /nc^) 

The relation between the work thickness and the optimum 
clearance was shown to be given by 


. 1+1. 5 e^ 

_ (^_.__ ) (2.+e^)+0.5e^+0.43 

O I 


( 1 . 6 ) 
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Fig.L9 Geometry of deformation of two 
fibres at a distance 8 apart 
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The penetration of the punch (A + CQ)/t was expressed as 


A c 

A + o 


0.43 + 0.5 


(1+1.5 e^)(l+e.) 

0.65+1.32 e^ ®f 


(1.7) 


The value of natural fracture strain = ln(e^+l) was 

substituted and the final expressions were given as 


t_ 

c. 


2.3 exp( ©£-1) 

= 1.36 exp( ep [ s - e^ ep r-f ^ 


( 1 . 8 ) 


and 


A + c 


T 1.9 exp( e^)-l 

— = — j 

^ 2.56 exp(C£)-l 


(1.9) 


The assumed pattern of the variation of punching force with 
punch travel is shown in Fig. 1.10. It is mentioned that the 
maximum punching force is attained just when the fracture starts 
and drops very rapidly with the progress of the punch. The 
maximum force (Fig. 1.11) was given by the following equation ; 


P 


max 



where L = punch width 

= average stress 

3 V 


( 1 . 10 ) 


From the same figure the work required, W, for the punching 
operation was expressed as 
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W 



max 


tp 


( 1 . 11 ) 


where 


(A + c^) 

p = the penetration r 


Then the values of the optimum clearance and penetra- 
tion depth (A + c^) in terms of the sheet thickness, for 
different values of ductility ( e^) of the material were cal- 
culated using the above equations. Some of the calculated 
values are given below : 


®f 

1 

1.5 

2.0 

2.5 

cjt 

0.215 

0.133 

0.082 

0.05 

(A + c ) 





t 

0.278 

0.289 

0.296 

0.30 


It was concluded from the results that depending on the 
ductility of the material the clearance varies from about 5 to 
20/ of the sheet thickness. It is suggested that for the blank 
ing of more ductile material optimum clearance required is less 
which was also observed by other workers in this field. Regard- 
ing the percentage penetration at fracture it was about 30^ of 
sheet thickness and increased very slowly with the ductility of 
the metal. For the verification of the proposed model, Popat 
[25] performed a large number of experiments using a single 
solid punch and variable clearance die set up. Two die halves 
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were mounted on two adjustable members running in a dovetail 
slide, the required clearance was set using the screw and nut 
assembly and feeler gauges. The punch force-punch travel 
pattern was obtained on the storage oscilloscope screen using 
a load cell below the die and a LVDT fixed onto the press 
ram. Further it is mentioned that to decide the value of 
optimum clearance, all the fractured pieces obtained at various 
clearances were examined simultaneously, the optimum clearance 
being indicated by the cleanness of the edge produced. A 
two-sided cutting was performed on sheet materials with 
different ductilities and thicknesses, at various clearances. 

It is mentioned that the following observations were made from 
the results of the experiments. 

(i) The optimum percentage clearance remains almost constant 
for any one material, since the fracture strain is independent 
of the thickness, and it changes only with the ductility of 
the material. The Fig. 1.12 shows the variation of optimum 
clearance with fracture strain, for the values obtained from 
experimental results and those obtained using equation (1.8). 
The nature of the variation for the experimental and theoreti- 
cal curves seems to match well. In the similar way the values 
of the percentage penetration obtained from equation (1.9) and 
those from the experimental results are plotted on the same 
graph, they are also found to be in good agreement. But while 
comparing the values of the maximum punch force measured 
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Fig. 1.12 Variation of the optimum clearence and penetration 
with fracture strain. 
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experimentally and calculated using equation (1*10), it was 
found that theoretical predicted values were higher than the 
experimentally found values. This was attributed to the fact 
that rupture of fibres is not simultaneous, which was an assum- 
ption while deriving the theory, but the rupture is continued 
during the finite interval of time. Close agreement between 
the theoretically predicted and experimental values was shown, 
after inserting some corrections using energy criterion. 

It has been noticed that from sixties onwards with the 
advent of High Energy Rate Forming (H.E.R.F.) methods, which are 
suitable for ambient and high temperatures, most of the research 
on blanking is diverted to fast blanking. The H.E.R.F. methods 
have offered the possibility of fast blanking with reduction of 
unit production times, improvements in the sheared edge of the 
blank produced and reduction in the amount of energy required 
to effect the plastic deformation during the process. References 
26 to 31 are devoted to the investigations on H.E.R.F. methods 
applied to blanking operation. 

Other works [32,33] have been devoted to manufacturing 
and other aspects. 

1.4 OBJECTIVE AND SCOPE OF THE PRESENT WORK 

The main objective of this work is to develop a finite 
element model for analysing the elastic-plastic deformation 
during blanking operation. The finite element method is selected 



33 


for the analysis because it is a versatile method and it is 
possible to predict with sufficient accuracy the distribu- 
tion of stress, strain, strain-rate, etc. for a conplicated 
problem. It is intended to prepare a displacement based finite 
element programme employing 'initial stress method' for this 
nonlinear problem. 

The punch-die clearance is the most important parameter 
of blanking process. When the clearance is adequate (optimum), 
the cracks initiating from the punch and the die corners will 
meet resulting in clean edges. The amount of the work required 
to be done for complete separation of the blank will also be 
minimum. In such processes there is severe plastic deforma- 
tion before fracture. Thus a suitable crack initiation cri- 
terion based on the value of a combination of the strains is 
to bo identified. Then a criterion for determining crack 
propagation direction will also be postulated which will be 
necessary for estimating the optimum punch-die clearance. 

In addition to the fractional punch— die clearance, the 
material properties and the sheet-rthickness are also considered 
to be important process parameters. During the analysis a 
parametric study is being undertaken to determine the effects 
of each one of the above on the maximum punch force required, 
punch penetration at crack initiation* optimum clearance etc. 
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To provide an overall idea about the deformation process 
the distributions of principal and effective strains in the 
clearance zone will be shown, with the progress of the punch. 
The deformation pattern and its variation with punch penetra- 
tion will also be presented. 

In order to verify the correctness and accuracy of the 
results obtained from the analysis, these are to be compared 
with the results of the previous experimental and theoretical 
works. A more extensive study is required to actually deter- 
mine the complete crack path with the progress of the punch 
and is not included in the present work. It will be shown 
that an attenpt is being made for getting the crack path. But 
it is felt that to obtain the correct crack path it requires 
different type of the element and a very fine mesh to dis- 
cretize the continuum. This will result in a very big size 
of the stiffness matrix and a large number of equations to be 
solved at each iteration. 

Although stresses, strains etc. are calculated assuming 
small deformation theory to be applicable, it is possible to 
take the large deformations into account by modifying the nodal 
coordinates, calculation of new stiffness matrix etc. after 
each displacement increment. Similarly with the necessary al- 
terations in the finite element model, the effects of the 
punch shear, frictional conditions at the tool-work interface. 
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clamping forces» strain rate etc, could be inplemented into 
the analysis. 

The present work is restricted to the two dimensional 
analysis of blanking operation. But in order to take into 
account the process of close contour blanking, a more rigorous 
three dimensional analysis should be performed using this 
technique . 





CHAPTER 2 


GENERAL THEORY OF FINITE ELEMENT METHOD 

2.1 INTRODUCTION 

The finite element method is much more versatile compared 
to other solution methods described earlier, and it is possible 
to analyse a complicated problem using this technique. There- 
fore, in recent years it has been also applied in the field of 
metal forming and sheet metal cutting [34 to 45] . With its use 
stress, strain, strain-rate and temperature distribution can be 
predicted more accurately, and hardness, ductility and grain 
homogenity can be deduced. 

This chapter briefly outlines the theory of finite element 
displacement formulation, as applied to the solution of problem 
of statics. Displacement formulation of elements used in the 
analysis of plane strain situation is described, using eight 
noded isoparametric elements. The two dimensional finite ele- 
ment displacement formulation of an elastic continuum is ex- 
plained. Methods of solution of plasticity problems using 
elasto-plastic analysis is stated in short. The two most 
commonly used yield criteria are described. The expressions for 
the elasto-plastic stress-strain matrix using von-Mises yield 
criterion are stated for the 3-dimensional as well as for plane 
strain situations. At the end, the steps during a typical 
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load/displacement increment is summarised for the elasto- 
plastic analysis using 'initial stress method'. 

2.2 FINITE ELEMENT METHOD 


In many situations an adequate engineering model of the 
system is obtained by replacing the system by a finite number 
of well defined conponents. Such problems are called dis- 
crete problems. In contrast to this are the continuous sys- 
tems in which the formulation leads to differential equations 
or equivalent statements, which imply an infinite number of 
elements [46]. Even if the number of elements is very large, 
the discrete problems can be solved readily now with the 
advent of digital computers. The various discretisation 
methods suggested from time to time to solve the realistic 
continuum problem both by engineers and mathematicians, all 
Involve approximations which approach the true continuum solu- 
tion as the number of discrete variables increase. 


The discretization of continuum problems has been 
approached differently by mathematician and engineers. The 
first have developed general techniques applicable directly 
to differential equations governing the problem, such as a 
^itir^ difference technique. On the other hand, engineers 
often approach the problem more intuitively by creating an 
analogy between real discrete elements and finite portion of a 
continuum domain. 
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It is for the engineering ’direct analogy* view that the 
term ’finite element’ has been born. Clough appears to be 
first to use this term, which implies in it a direct use of 
standard methodology applicable to discrete systems. 

Finite element method (FEM) is now very well described in 
a number of books [46 to 50] . However, it is briefly outlined 
here to keep continuity of the discussion. FEM is a numerical 
discretisation procedure by the use of which a wide range of 
complex boundary value problems can be analysed. Elastic, non- 
linear elastic and elasto-plastic constitutive relations can be 
implemented within finite element framework in a straight for- 
ward manner. 

In FEM, the region of interest is divided into a finite 
number of simply connected sub-domains or elements. An approxi- 
mate functional value of the solution is assumed over the sub- 
domain so that the parameters, say u^^ (i = 1, n) of the 

function becomes the unknown of the problem. 

2.2.1 Physical Interpretation of FEM 

The basic concept is derived from structural analysis. 

Here every structure is approximated as a physical assemblage 
of individual structural components or finite elements. The 
elements are interconnected at a finite number of nodes and 
sometimes along the boundaries of the element. Assuming the 
■approximate behaviour of individual elements, the behaviour of 
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the entire system can be analysed. After assembling the 
individual elements, the necessary boundary conditions are 
imposed, with primary nodal values as unknowns* Solution of 
the resulting set of equations yield the system response. 

Considering the nature of each individual component in 
an element, the relationship between the primary variable u§^^ 
and the forcing function Q^, for a linear system can be 
written as [46] 


Q? = Kf . uj 4- Q^, 


( 2 . 1 ) 


repetitive indices indicating summation. 

Connection among elements can be established using, 


(a) one set of variables, u ^ , for the assembled system, i.e., 
the condition of continuity 


u- = u? (2.2) 

kJ sJ 

( b) the equilibrium of the second set of variable Q®, at 
each node and equating it to zero, 

NEL ^ 

Z Qf = 0 (2,3) 

e=l 

in which, NEL, is the number of elements considered. Using 
equations (2.1) to (2,3), the resulting equations can be 


written as 
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[k] £u ? + £Q} = 0 


in which 


NEL 


[K] = K, . = Kj. 


and 


NEL 

fQ] = = E qJ3_ 

e=l 


(2.4) 

(2.5) 

( 2 . 6 ) 


2.2.2 Finite Element Displacement Formulation of an Elastic 
Continuum 

The special features of the problem are given below [51] 

( a) The continuum is separated by imaginary lines or surfaces 
in a number of 'finite elements'. 

(b) The elements are assumed to be interconnected at a dis- 
crete number of nodal points situated on their bounda- 
ries. The nodal displacements will be the basic unknowns 

(c) A set of functions is chosen to define uniquely the 
state of displacement within each finite element in terms 
of nodal displacements. 

(d) The displacement functions now define uniquely the state 
of strain within each element in terms of nodal displace- 
ments. These strains, together with any initial strains 
and constitutive properties of the material will define 
the state of stress throughout the element and hence 
along the boundaries . 



(e) A system of forces concentrated at the nodes and equilibra- 
ting the boundary stresses and any distributed load, is 
determined. 


Once this stage is reached, the solution procedure follows 
the standard discrete system pattern. The above procedure 
involves a series of approximations. Ensuring displacement 
compatibility between adjacent elements may not be possible 
always, though within each element it is obviously satisfied. 
Concentrating the equivalent forces at nodes, equilibrium 
conditions are satisfied in the overall sense only. 


The displacements can be expressed as 

e 


u 


^i| 


E N^a? = .,.,3 = N a 


(2.7) 




where u = displacement at any point within the element 

N = prescribed function of position called shape function 
a® = nodal displacement vector for a particular element. 

The strain vector, e , can be obtained with displacements 
known at all points within the element. Dropping subscript for 
element, 

e = B a (2.8) 


For a linearly elastic constitutive law. 
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where 

= initial strain matrix 
D = symmetric matrix 


— initial residual stresses. 


q 


e 



= nodal force vector which are 
statically equivalent to the 
boundary stresses and distributed 


loads on the elements . 


To make the nodal forces statically equivalent to the actual 
boundary stresses, and distributed loads, the simplest procedure 
is to impose an arbitrary virtual displacement and to equate the 
external and internal work done by the various forces and 
stresses during that displacement (i.e.) 

6a® q® = 6a® [ / ad(vol) - / N^bd(vol) 

V® V® 

- / t d(Area)3 (2.10) 

A® 


where 

b = vector of body force components 
t = vector of traction components. 



43 


As this relation is valid for any value of virtual displacement. 


q® = / ad(vol) - / N^bd(vol) - / t d(Area) (2.11) 


V'= 


V" 


qG = K®a® + f® 


(2.12) 


K® = / B^DBd(vol) 

,re 


(2.13) 


f® = -/ bd(vol) - / t d(Area) ~ / B^ D e d(vol) + 

V® A® V® 


V 


f d(vol) 

e 


(2.14) 


External concentrated forces may exist at the nodes 
denoted by vector 


‘2 


I 


(2.15) 


Equation (2.7) when applied to the whole structure will 


give] 


u = N a 


(2.16) 


in which a is nodal displacement vector for all the nodal points 



44 


Similarly, 

e = B a (2.17) 

Dropping the bar superscript 

6a^r = / 6 u^b dV + f 6 J t dA - / 6 cj dV (2.18) 

V A V 

Substituting, equations (2,16), (2.8) and (2.9), once 
again a system of linear equations are obtained. 

Ka + f = r (2.19) 

K = / DB dV (2.20) 

V 

f = - / N^b dV - / t dA - / B^D s. dV + / B^ a dV (2.21) 
V A V ° V ° 

The integrals are taken over the whole volume V and over the 
whole area A on which the tractions aro given. 


K. j 

fi = 2 fj 

(2.22) 

/ N’^b dV = 2 

S N^b dV 

(2.23) 

V 



/ t dA = 2 

/ t dA 

(2.24) 

V 

A^ 



The same is true for other elements in the equation (2.20) 
and (2.21), 
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2.3 DISPLACEMENT FORMULATION OF THE ELEMENT USED 

The unique description of the displacement within each 
element in terms of nodal values at boundary points or internal 
points of the element is the basic step in any displacement 
finite element formulation and can be expressed as 

u = N a® 

where N is the matrix of shape functions. 

If the shape function chosen to describe the element 
geometry are identical to those used to prescribe function 
(displacement) variation, then the element is termed iso- 
parametric. The basic procedure in the isoparametric finite 
element formulation is to express the element coordinates and 
element displacements in the form of interpolation functions 
using natural coordinates system of the element. The fun- 
ctional property of the interpolation function N^^ is that its 
value in the natural coordinate system is unity at node i and 
zero at all other nodes. The formulation given below follows 
general pattern of derivation suggested by Bathe and Wilson 
[47] and Zienkiewicz [46] . For a two dimensional element 
interpolation functions for eight noded case is given below. 
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N;L =1 (l-s)(l-t^) - (i) N5 - (|) Ng 

N2 = I (l+s)(l-ti) - (|) • (|) Ng 

Ng = I (H-s)(l+t^) - ( 1 ) - ( 1 ) 

N4 = I (1- s)(l+tj_) - (|) - (|) Nq 

% = I (l-s2)(i_t^) 

Ng = ^ ( 1— tj^) ( 1+ s) 

= i (l--s2)(i+t^) 

Nq = I (l-t^)(l-s) 


(2.25) 


A quadratic variation of displacement and geometry are 
assumed for this element and it has 16 displacement degrees 
of freedom (Fig. 2.1). 

Cartesian coordinates x,y and natural coordinates s, t^^ 
are related as follows ; 

X = Nj^x^ + N2X2 + NgXg + N^x^ + ^6^6 ^.^^7 "*■ ^3^8 

(2.26a) 

y = Nj^y^ + N2y2 + Ngyg + N^y^ + N^y^ + N^y^ + N^y^ + Ngyg 

(2.26b) 

where Nj_ to Ng are given in equation (2.25) . 
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(a) Local co-ordinate system (b) Natural co-ordinate 

system 



4 fx4 

(c) Parabolic variation of surface traction on 
straight edge 43 

Fig. 2. 1 Eight noded parabolic isoparametric element 
for two dimensional analysis 
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Similarly, the displacements expressed in terms of nodal 
values are ; 


u + N^u^ + 

V = + N2V2 + N3V3 


+ N4U4 + + NgUg 


+ N4V4 + N5V5 + NgVg + N7V7 + NgVg 


(2.27a) 

(2.27b) 


Along the edge of the element the variation of displacement is 
quadratic . 

VVriting equation (2.27) in matrix notation 
{u} = [N] {a} 

(a ! T = UjVj U3V3 UgVg] ( 2 . 28 ) 

Element strain vector for plane strain case is 

5®’ ^ = C "xx (2-29=>) 

whore Exx = 37 ’ ®yy " ^ ’ ^xy = 37 + 3^ (2.29b) 

The strains at any point within the element is obtained 
knowing the nodal displacement, by 

{e} = [b] £a} (2.30) 


where B is the strain displacement matrix given by 
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[B] = 


dN, 

0 

dN2 


dNg 

0 

dN4 

0 

dNg 

0 

dNa 

0 

dN„ 


dNg 


W 

dx 

0 

dx 

dx 

dx 

dx 

dx 

0 

dx 

0 

0 

dN. 

- 0 

SNj 

dy 

0 

dNg 

dy 

0 

SN4 

dy 

0 

dN^ 

0 

dN^ 

dy 

0 

dN^ 

dy 

0 

SNg 


dNjL 

dN2 

dN2 

dNa 

dNg 

dN^ 

dN4 

dNs 

dNj 

dNg 

dNg 

dN^ 

dN^ 

dNg 

dNg 

dy 


dy 


W 


dy 

dx 

dy 

dx 

dy 

dx 


dx 

dy 

dx _ 


(2.31) 


For plane strain the elasticity (stress-strain) matrix can be 
written as 

*~1-Y Y 0 


[D] = 


(1+y)(1-2y) 


Y 

0 


1-Y 

0 


0 

1-2y 

2 


(2.32) 


The matrices to be evaluated are the stiffness matrix K, the 
body force vector, and the surface load vector for distributed 
loads at the boundary. The matrices in general are evaluated by 
numerical integration, using Gaussian quadrature formula. For 
the plane strain condition thickness being unity, the matrix K 
can be obtained as 

n n 

[K] = t S 2 a. a F. . (2.33a) 



T 

b: . 

13 


D B. . 
13 


det J. . 


(2.33b) 
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in which n is the number of integration points in each 
direction, det J. . is the determinant of the Jacobian matrix 

-J 

of transformation 

^i » *^3 Gaussian quadrature weights 

Fj^j are the function values to be evaluated at sampling 
points s^, t^j , 

For the evaluation of stiffness matrix of the element, a 
2x2 numerical integration is used. The function values are 
evaluated at 4 sampling points and summed upto get the stiffness 
matrix of the eight nodes element* 

2.4 ELASTO-PLASTIC ANALYSIS 

The difference in elastic and plastic behaviour under 
uniaxial stress is clearly brought out in Fig. 2.2. In non- 
linear elastic behaviour, the stress can be expressed as a 
function of strain as o = a(- e ) . The main difference of 
plasticity formulation from nonlinear elastic formulation is, 
such an explicit relationship is not evailable. Although the 
stresses at any level of strain have to lie on or within the 
current yield surface, the exact value of each component 
cannot be determined. 

2.4.1 Development of Numerical Techniques for Plasticity 
Analysis 

Theoretical as well as practical interest in plasticity 



a. Linear elastic b. Non-linear elastic 


a 
a 


c. strain hardening d. Ideal plasticity 
plasticity 

Fig. 2.2 Elastic and plastic behaviour under 
uniaxial stress 
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has stimulated much work in the numerical solution of 
elasto-plastic problems by the FEM and other procedures [52] . 
The literature on small strain displacement elastic-plastic 
solution using FEM is now large [53 to 63] . Only the 
difference lie between the apparent forms of the nonlinear 
solution algorithms used, and the form of constitutive rela- 
tions postulated. 

The common conputational approach used with incremental 
theory of plasticity in finite element method are t 

(i) incremental tangent stiffness approach, 

(ii) 'initial strain' method, and 

(iii) 'initial stress* method [48,62]. 

In the first method the change of material constitutive matrix 
forming element stiffness, assembling and inverting an overall 
stiffness matrix for each step of analysis is to bo carried 
out. In 'initial strain' and the 'initial stress* approaches, 
in all iterations a simple resolution of constant linear 
elastic problem is carried out. The 'initial stress' method 
is suitable for solution of ideal plasticity or for strain 
hardening formulations. 

2.4.2 Yield Criteria 

A yield criteria is a hypothesis concerning the limit of 
elasticity under any possible combination of stresses. The 
suitability of any proposed yield criteria must be checked by 
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conducting experiments. 

If a point in a ductile material is subjected to the 
principal stresses , cr^, and is represented by a Mohr 
Circle, and if the principal stresses are changed to 
( and then the size of Mohr circle, for this 

new state of stress, remains same but is shifted by a distance 
along the e-axis. The additional stresses make up a 
hydrostatic (tensile or compressive) stress system. It is 
found that the absolute size of the Mohr circle alone deter- 
mines the limit of elasticity and is independent of its 
position. This is to say that the yield criteria is a fun- 
ction of and independent of 

hydrostatic stress component (cj^+cr2+^3)/3. , 

Thus yielding occurs when some scal^ function of the 
principal stress differences reaches a critical magnitude, 
mathematically : 

~ constant (2.34) 

Tresca Yield Criteria 

Perhaps the simplest function imaginable which satisfies 
equation (2.34) is of the form 


O'. 

1 


Oj = constant 


(2.35) 
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Equation (2,35) interprets when largest of the three 
magnitudes ® critical constant 

value (for a given material) yielding begins. 

This was first suggested by Tresca in 1864. It states: 
'yielding occurs when the greatest absolute value of any one 
of the three maximum shear stresses in the material reaches a 
certain value* . 

Applying this criteria for a state of pure shear and 
simple tension following relation could be derived : 




= o 


yp 


= 2k 


(2.36) 


where o is tensile yield stress and k is a yield shear stress, 
y H 

Von-Mises Yield Criteria 

Another admissible function, which satisfies equation 
(2.34) is of the form 

2 2 2 

(Oi~a2) + (a^-o^) + (o^-o^) = constant (2.37) 

In this type of function, each of the principal stresses con- 
tributes to yielding. This function was proposed by Huber 
(1904) and von~Mises (1913) and by J.C. Maxwell in a letter to 
Kelvin in 1856 [64], The criteria was then interpreted by 
Hencky as 'yielding begins when the shear strain energy 
reached a critical value' . 
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The value of the constant in equation (2*37) can be again 
determined from the simple tension yielding, i.e., = cTyp^ 

02 = 0, o^ =0 and pure shear yielding, i.e., o^ = -o^ =k , 

O 2 = 0. Hence yield criteria becomes 

+ (o^-o^)^ = 2o^ = 6k^ (2.38) 


For a general state of stress (o^, o^, ^z»‘’'xy» "^yz’ *’'zx^ ’ 

the von-Mises criteria becomes 


[( 0-0 )^ + (o- 

J2 X Y y 


o )■ 
z 


+ (a_ 


■a )■ 

X 


+ 6 ( 

^ xy yz 


)f= 

zx"^ ■’ 


o 


yp 


(2.39) 


In Fig. 2.3a and Fig. 2.3b both the yield criteria are represen- 
ted in principal stress space and tc plane respectively. 


2.4*3 Elasto-Plastic Stress-Strain Matrix 


Hooke’ s law for the isotropic elastic material can be 
written in matrix form as 


ia] = [D] £e} (2.40) 

where [d] is given by equation (2.32) for the plane strain 
condition. 

By now it has been well established that the plastic defor- 
mations are incremental. Yamada and Yoshimura [53] have deve- 
loped a plastic stress-strain matrix (Elasto-Plastic Matrix) 
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by inverting the Prandtl-Ruess equations obeying the von-Mises 
yield criterion. This symmetric plastic stress-strain matrix 
is given by following equation 


D®P] = 



(2.41) 


lere 


2 7^ H’ X - 

30 ( 1 + 3q ) , 0 


’k + (0 -0 + (0 -cr^)^ 

2 ^ X y'^ '' y '■ z x'^ 


+ ■'xy 1 


or the plain strain case, the [D®P] will reduce down to the follow 
ag matrix : 


CM Nja) 
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izl 

1-2y 



S 


[D®P] 


E 

1+Y 


Ji 

1-2y 


a ’ a ' 
X y 







(2.42) 


where 

S = I -02 (1. ^) = § -^(1. . 

2.4.4 The 'Initial Stress' Computational Method 

In the 'Initial stress method' [5l] the solution of the non- 
linear problem is approached in a series of approximations. During 
a load/displacement increment, a purely elastic problem is solved, 
determining an increment of strain lA^ and stress 

at every point of the continuum. For the increment of strain 
found from the analysis, the stress increment in general will not 
be correct. If iAa} is the true increment of stress possible 
for the given strain tl^^ the situation can only be maintained by a 
set of body forces, equilibrating the initial stress system given 
by {Act' 3 - {AoJ . At the second stage of the conputation, 

this body force system can be removed by allowing the structure, 
(with unchanged elastic properties) to deform further. An addi- 
tional set of strain and corresponding stress increments are 
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caused. Once again these are likely to exceed those permi- 
ssible by the nonlinear relationship and the redistribution of 
equilibrating body forces has to be repeated. If the process 
converges then within an increment the full nonlinear conpati- 
bility and equilibrium conditions will be satisfied. In each 
cycle, since the same elastic problem is solved only partial 
inversion of the elastic equations is needed. 

For the elastic-plastic situation, the steps during a 
typical load/displacement increment can be summarised as 
follows [ 56 ] : 

!• Increment of displacement is applied. Elastic increments 
of strains {Ae‘}^ and elastic increments of stress 
which corresponds are found. 

2 . is added to stresses existing at start of 

increment, I^q} » to obtain £0'} , and thus o* , Check 

whether 0* < If the condition is satisfied only 

/ -H 

elastic strain changes occur and the process is stopped. 

If the condition is not satisfied proceed to the next 
step. 

3 . If a* _> 0^ and if 0^ ~ ( the element has yielded 
at start of iteration), £^0}^^ is found by 

£A 0 1 1 = £A e» where [D®^] is computed using 

I0’ } . Stresses which are supported by body forces 


are evaluated. 



60 


= £Ac'3^ - CAa3j^ 


(2.43) 


Current stress is 
£a} = Cu’ }+£ 


(2.44) 


Current strain is 


Is} = {e’ } + { Ae’jj^ 


(2.45) 


4. If a* > CT , but cr^ < a , i.e., elastic-plastic 

r65l 

transition is occurred in that iteration, a factor L J 
equation (2.46) is found and then the interniediai® stres 
value is found out at which yield begins. The increment m 
stress is computed by the relation, £ Ao [D ^3 ® ^1 

starting from that point. Then proceed as in step 



Factor = 



(2.46) 


* T *1 "t’ 

Nodal forces are computed corresponding to the equi 
ing body forces. These are given by 


{ Rl^I = cl(vol) 


(2.47) 


Elemental nodal force are added to get the complex® nodal 
force vector . 
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6. The problem is resolved using the original elastic pro- 
perties and the load system { to find {Act*} 2 

{A s*}2. 

7. Current values are found. 

8. Steps 2 to 6 are repeated. 

The cycling can be stopped when the appropriate conver- 
gence criterion specified is satisfied. 

2.4.5 Convergence in Finite Element Method 

Since finite element method is an approximate numerical 
procedure for solution of engineering problems, like any other 
approximate numerical techniques, the question of convergence, 
rate of convergence and accuracy of the method are of great 
importance. 

During the procedure in developing finite element method 
employing variational, residual or virtual work methods, the 
key equations involve integrals over the domain (and also in 
some cases over the boundary), it is seen from equations (2.20) 
and (2.21). If each of these integrals tends to its exact 
value as the element size decreases to zero the finite element 
method would yield solution \Miich correspondingly tends to the 
exact solution. So the requirements for any displacement for- 
mulation can be stated as follows [48] : 
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( a) The displacement model must be continuous within the 

elements and the displacements must be compatible between 
the adjacent elements. 

(b) The displacement model should include the rigid body dis- 
placements of the element. 

(c) The displacement model should include constant strain 
states of the element. 

Element formulations which satisfy criterion (a) is called 
compatible or conforming and formulations which satisfy (b) and 
(c). are called complete. 

The adoption of an incremental plasticity theory and itera- 
tive solution technique, such as used in 'initial stress 
method' loads to errors due to size of increment and finite 
values of criteria for convergence specified. Convergence 
criteria specified in a solution, indicate whether the non- 
linear equations describing the behaviour of the material when 
subjected to increment of load/displacement have been solved 
with sufficient accuracy. This is achieved by comparing the 
change in some aspect of the solution between iterations and 
terminating the procedure once this change has become suffi- 
ciently small. 

Nayak and Zienkiewicz [6l], and Rowe and Davis [66] 
have described some of the convergence criteria for load based 



analysis and displacement based analysis. The convergence 
criteria imposed in this work is on the change in the absolute 
magnitude of largest term in load correction vector between 
the iterations, compared as a percentage of the absolute 
magnitude of largest term in load correction vector at that 
instant calculated by equation (2.47) . 



CHAPTER 3 


FINITE ELEMENT ANALYSIS OF BLANKING 

3.1 INTRODUCTION 

This chapter describes the finite element formulation of 
sheet blanking operation, which is considered as a two dimen- 
sional process occurring under plane strain conditions. Eight- 
noded isoparametric quadrilateral elements are used for the 
analysis and a 2x2 numerical integration is employed for cal- 
culation of stresses and strains . A list of the assumptions 
made and the description of the data vhich is necessary for the 
execution of the programme is also mentioned, together with the 
solution procedure. The criteria used to determine crack ini- 
tiation and crack propagation direction are also stated. The 
finite element code was verified by a punch indentation pro- 
blem [61]. A short description of the main programme and 
important subroutines is given in Appendix A. 

3.2 PROBLEM DESCRIPTION 

The actual blanking operation is normally a closed con- 
tour sheet metal cutting process involving a 3-dimonsional 
deformations which is rather difficult to analyze. For the 
present analysis the process is simplified as a two dimensional 
under plane strain conditions, just similar to the two sided 
bar cropping operation. In a simple blanking operation 
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(Figure 3.1) a strip of thickness t is cut along the sides AA 
and BB with a sharp edge punch and die set having the clearance 
c. As the punch moves down the sheet metal deforms, initially 
the deformation being elastic in nature and with the further 
movement of the punch, material deforms plastically. When the 
penetration of the punch reaches a certain critical value, 
cracks start at the die and the punch corners. With the fur- 
ther movement of the punch the cracks propagate inside even- 
tually leading to the separation of portion of the sheet metal 
under the punch (knowi as a blank), from the rest of the strip. 
It is known that the severest deformations occur in the 
clearance gap. In the present work an attempt has been made 
to obtain the distribution of stress and strain fields in the 
clearance gap so as to determine the crack path and the varia- 
tion of punch load with punch travel. 

3.3 FINITE ELEMENT FORjVlULATION 

The process of sheet blanking is analyzed using the 
finite element method with the following assumptions. 

(i) The blanking process is a two dimensional occurring under 
the plane strain conditions. 

(ii) The process is quasi-static with the punch moving at slow 
speed and hence the effect of the strain rate is neglected. : 




Fig,3.l Schemetic diagram of blanking operation 
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(iii) The bottom surface of the punch and the top surface of 
the die are flat and smooth, and their' edges are sharp. 

(iv) The cutting conditions are identical at the punch and the 
die edges. 

(v) Small deformation theory is applicable for the purpose of 
calculating stresses and strains. 

(vi) The material of the sheet is isotropic, linear elastic 
linear strain hardening type obeying von-Mises yield criterion. 

Fig. 3.2 shows the two dimensional blanking process. The defor- 
mation pattern is symmetric about the centre line X-X and, 
hence, half of the sheet needs to be analysed, say region ABCD 
(beyond the deformation in the sheet is assumed negligible). 

The detailed model of the region ABCD is shown in Fig. 3.3, in 
which for convenience two rigid blocks separated by the dis- 
tance c (clearance gap), are moving in opposite directions. In 
the present work displacement controlled analysis is performed 
and at each increment of displacement the two blocks move equally 
in opposite directions. 

The following geometric boundary conditions are specified 
for the present work (refer Fig, 3.3). 

(i) The boundary sides AB and CD are fixed in the x-direction 
and free in the y-direction. 




Punch 
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(ii) The boundary side EA is free in the x-direction and uni- 
form downward displacement is specified in the y-rdirection. 

(iii) The boundary side CG is free in the x-direction and uni- 
form upward displacement is specified in the y-direction. 

3.3.1 Element Type and Step Size 

A general purpose elastic-plastic finite element programme 
using ’initial stress method' was developed on DEC-1090 computer 
system for the analysis of the present problem. 

The element selected for the discretization of the continuum 
was eight noded two dimensional isoparametric element which is 
already described in Sec. 2 . 3 . ■ For this type of parabolic 
element the variation of displacement is quadratic within the 
element and hence the variation of stresses and strains is 
linear. 

In this incremental cum iterative technique of elastic- 
plastic finite element analysis the size of the increment and 
the number of iterations to be performed within an increment 
will affect the results obtained. For this work the values of 
the displacement specified in first displacement step is such 
that yielding just begins. Because of the antisymmetry, two 
gauss points yields whoso effective stress just equals the 
yield stress of the material. It has been observed [46] that 
by using the 'initial stress method' for the nonlinear material 
problems, even if the large values of load/displacement 
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increments are specified, the results obtained are quite satis- 
factory. From the kind of the mesh selected for this analysis 
the total number of degrees of freedom comes out to be around 
800, which requires large computer storage and the longer exe- 
cution time for the programme. Thus it was decided to give 
the total displacements required to be applied (till the indi- 
cation of the crack initiation is obtained), in about 5 to 8 
steps, depending upon the local fracture strain value specified. 

3.3.2 Criteria for Crack Initiation and Crack Propagation 
Direction 

Shaw and Avery [24] have mentioned that metal forming and 
sheet cutting (i.e., cold working) operations involve severe 
plastic deformation in general, i.e., the material is strained 
much above the elastic point. The sheet blanking operation is 
also one of the cold working processes, during which it is 
observed that with the penetration of the punch into the sheet, 
the metal lying between the punch and die (in the clearance gap) 
deforms severely and finally the crack develops and propagates. 
Crack initiation is a local phenomenon and here it is postula- 
ted that crack is initiated at any point where the effective 
strain ( © ) just exceeds the local fracture strain value ( e^) 
of the material. In the present analysis the effective strain 
is calculated, from the values of the individual strains, using 
following, expression for 2-D plane strain case 
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r= 2/3 [ + Syy + 0.5 (3.1) 

This & was then compared with the local fracture strain, 
of the material, i.e., if ^ >. then the crack is assumed 
to be initiated at that point. If e < at all the points 

then the next increment of displacement is applied and process 
repeated. 

As the shearing strains are predominant in the blanking 
process, blanking is also referred to as a process of shearing. 
This was observed from the results (presented later in 
Chapter 4) of the finite element analysis, i.e., the shearing 
strains comes out to be higher than values of the normal 
strains. Thus here, it is postulated that the direction of 
the crack propagation is the same as the direction of the 
maximum shear strain at crack initiation point. 

3.4 INPUT DATA STRUCTURE 

The required finite element mesh was initially generated 
using an automatic mesh generation programme which is a revi- 
sed version of the Durocher and Gasper's [67] programme. The 
following data was supplied as input to the eiasto-plastic 
programme for the purpose of analysis. 

(i) Total number of nodes and elements, number of nodes per 
element, number of displacement prescribed boundary nodes, 
order of numerical integration, number of degrees of freedom at 
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each node, maximum number of iterations for each displacement 
step, bandwidth of the matrix etc. 

(ii) Material properties, i.e., the values of the elastic 
modulus, Poisson’s ratio, yield stress, strain hardening co- 
efficient and local fracture strain. These values are taken 
from references [25] , and [34]. 

(iii) Nodal data which includes the node number and, x and y 
coordinates of the node, for all the nodes. 

(iv) The element data in which the element number is given and 
the nodal connectivities of the element, for all the elements. 

(v) The boundary node number on which geometric boundary 
condition is known or prescribed and the type of the boundary 
condition, i.e., either free, fixed or some prescribed displace- 
ment, for all nodes having geometric boundary conditions. 

3.5 PROGRAMME TESTING 

The solution algorithm of the general purpose elasto- 
plastic finite element programme using ’initial stress method' 
is already described. The programme which was developed is 
applied initially for the solution of a punch rndentation pro 
blem [61], whoso results are known. Fig. 3.4 shows a rigid 
punch indenting into a metal specimen whose material properties 
have been varied from isotropic strain hardening (H’ = O.IE), 
ideally plastic (H' =0), to strain softening case{H’= -O.IE). 
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The problem being symmetric about the centre line X-X of the 
punch, only half the continuum ABCD is analysed. For this pur- 
pose the region ABCD is divided into 107 nodes and 28 parabolic 
isoparametric elements as shown in Fig. 3.5. For every incre- 
mental displacement step the mean punch pressure has been cal- 
culated, from which the pressure-displacement curves are drawn 
as shown in Fig. 3.6, for various values of H’ . As it can be 
seen from the graphs in the said figure, that the results of 
the present analysis match with those given by Nayak and 
Zienkiewicz [61] . For the case of ideally plastic situation 
the collapse load predicted by slip line solution [68] also 
matches with the finite element solution, which can be seen 
from the same figure. 



Fig. 3.4 Plane strain punch indentation problem 



Fig. 3.5 Discretization of region ABCD 
(of Fig. 3.4) 



0.001 0.002 0.003 0.004 0.005 0.006 0.007 

8/b 

Fig. 3.6 Punch pressure- punch displacement curves for 
punch indentation problem . 



CHAPTER 4 


RESULTS AND DISCUSSION 


4.1 INTRODUCTION 

This chapter presents an analysis of the results obtained 
for the blanking problem using the elasto-plastic ‘initial 
stress’ finite element programme. The continuum is discre*“ 
tized into 396 nodes and 117 eight noded quadrilateral ele- 
ments. The main variables for this problem are the material 
properties E>Y»^yp »H' and (the material is assumed to 
remain isotropic during the process), the sheet thickness t 
and the fractional punch-die clearance c/t. In the present 
work a large number of computer runs are taken , keeping the 
values of E = 12500 x 10*^ N/m^, y =0.34 and = 8 x lo"^ N/m^ 

as constant throughout the analysis. The values of other 
variables are changed i.e., H* is varied from 25 x 10 N/m to > 
100 X lo”^ N/m^, from 1.1 to 1.59, t from 1.47 x lO”^ m to 
6.0 X lO" m and c/t from 0.03 to 0.20 respectively; at a time 
only the value of one variable is changed. The CPU time re- 
quired for one run on the average comes to about 40 to 45 minu- 

i 

tes on DEC- 1090 computer system. Initially the effect of the ! 

various variables on the punch force-punch penetration pattern ! 
is analysed. Then from these results the fractional punch pene- ! 
tration ( — required at crack initiation is obtained for all | 
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the cases and these are discussed. The punch force required to 
initiate the crack is taken as the maximum punch force, 
and effects of various parameters on it are studied. From the 
values of the principal strains and their directions, for 
various c/t ratios and an indication of the optimum clea- 

rance ratio is obtained and the dependence of the optimum 
clearance ratio on the ductility of the material ( is 

brought out. Finally an attempt is made to explain the mode 
of metal deformation and distribution of strains, at various 
stages of punch penetration, in the punch-die clearance gap. 

The convergence criterion specified is, the change in absolute 
maximum value of residual load vector between iterations as one 
percent of its absolute maximum value, in the present work. 

4.2 RESULTS AND DISCUSSION 

Fig. 4.1 shows the finite eleinent discretization of the 
region ABCD of Fig. 3.3 together with the geometric boundary 
conditions. As mentioned earlier the discretization is obtain- 
ed using automatic mesh generation programme. The same mesh 
was employed for the conplete analysis of the problem but with 
the numerical values of thickness (t) and clearance (c) varied 
as and when required. This mesh is selected after checking the 
results of some sample runs with different meshes, A mesh finer 
than the used here was not selected, because that would lead to 
more number of degrees of freedom (with the present TOsh the 
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total degrees of freedom is about 800) . This means with finer 
mesh there is higher band width of the global stiffness matrix, 
larger number of equations to be solved at each iteration etc., 
thus requiring higher computation time. The size of the mesh 
is also restricted because it is decided to take a good number 
of computer runs by varying different parameters. 

Displacement controlled analysis is performed i.e., equal 
and opposite vertical displacements 6/2 are specified on sides 
AE and GC (Fig. 4.1) of the continuum, during each displacement 
increment step. The output of the computer run gave the value 
of the equivalent vertical load P/2 corresponding to 6/2. 

Thus, the value of load P (called punch force) was obtained 
against the percentage displacement ((6/t)xl00) also called 
the punch penetration, at each increment. The increments of 
applied displacement is continued till the crack initiated in 
any of the element (i.e., the crack initiation criterion 
stated earlier is satisfied). The punch force corresponding 
to percentage penetration upto crack initiation (( — )xlCX)) 
is called maximum punch force. From these results the varia- 
tion of punch force-punch penetration and other plots are 
drawn. 

4.2.1 Punch Force-Punch Penetration Characteristics 

Fig. 4.2 shows a set of such plots where the percentage 
clearance (| x ,100) is taken as 3,6,8.5,14 and 20, respectively. 
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with the values of sheet thickness, strain hardening parameter 
and local fracture strain taken as 1.47x10“^ m, 25x10”^ n/n? and 
1.25, respectively (Table 4.1a). Figures 4.3, 4.4 and 4.5 show 
similar sets of plots as the previous one with the values of 
as 50x10 N/m , 75x10 N/m and lOOxlo”^ N/m^ respectively 
(Tables 4.1b, 4.1c and 4*ld). In Figures 4,6 and 4.7 the value 

of thickness for a set of variation of percentage clearance is 

“3 -3 

taken as 3x10 m and 6x10 m respectively (Table 4.2), keeping 
the other parameters constant. 

From the above plots it is observed that within the range 
of clearances used, with the increase in percentage clearance, 
for the same amount of punch force the punch penetration in- 
creases. Again with the same amount of percentage clearance, 
as the H' is increased the curve moves upwards, i.e., the punch 
force required for same amount of penetration of the punch, 
increases. This is because with higher H' , the material is more 
strain hardened and its ductility reduces, thus, demanding more 
force for plastic deformation. The similar observation is made 
(i.e., the punch force increases) with the increase in sheet 
thickness. This is because the amount of material in the 
clearance gap increases as the thickness is increased and so 
requires more force. 

4.2.2 Punch Penetration at Crack Initiation 

In Fig, 4.8 a set of curves is shown in which the percentage 
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Table 4.1 Punch Penetration vs Punch Force at Different 
Fractional Clearances (t=l. 47x10-3 m) 

E = 12500x10^ n/m , y =0.34, = Sxio'^ N/m^, = 1,25 

a. H’ = 25x10^ N/m^ 


c/ t=0.03 

c/t=0 

.06 

c/ "fc-rO • 085 

c/t=0.14 

c/ t=0 .20 


{6/t)x P^- - 
100 ‘ 

(6/t)x 

100 

Pn 

(6/t)x 

100 

Pn 

(6/t)x P 
100 

(6/t)x P 

100 


0.02 33.3 

0.02 

33.0 

0.02 

33.8 

0.02 34.0 

0.03 34.3 


0.11 106.6 

0.12 

125.4 

0.13 

125.0 

0.15 123.2 

0.16 121.8' 


0.57 208.2 

0.63 

186.4 

0.67 

173.6 

0.75 159.0 

0.84 153.9 


2. '88 239.5 

- 3.18 

218.6 

3.39 

207.4 

3,78 188.6 

4.24 176.1 


14.45 335.0 

15.93 

317.4 

16.99 

308 .2 

18.91 288.4 

21.22 275.2 


23.67 354.0 

28.69 

364.8 

30.59 

353.6 

35.64 343.8 

38.25 320.6 



b. H’ = 50x10^ n/rt? 


c/t=0 .03 

c/t=0.06 

■“ c/t=0.085 

c/ 1=0.14 

c/t=0.2d“’ 

(6/t)x P 

100 

X 


(6/t)x 

100 


(«/t)x P 
100 

(6/t)=. 

100 


0.02 31.3 

0.02 

33.0 

0.02 

33.8 

0.02 34.1 

0.03 

34.3 

0,11 128.3 

0.12 

126.7 

0.13 

126.2 

0.14 124,4 

0.16 

102.8 

0.57 214.7 

0.63 

198.3 

0.67 

181.2 

0.75 166.4 

0.84 

160.2 

2.88 263.9 

3.18 

247.2 

3.39 

234.8 

3.77 216.2 

4.24 

203 .4 

14.45 457.4 

15..93 

448.0 

16.99 

434.2 

18.91 418.2 

21.25 

404.2 

23.34 505.8 

28.58 

531.2 

30.59 

544.6 

35.64 546.8 

41.82 

556.4 


contd 


• ♦ ♦ 


H’ = 15x10^ N/m^ 
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c. 


c/t=0.03 

c/t=0 

.06 

c/ t=0.085 

c/t=0. 

14 

O j 

1 

!1 

O 

.20 

6/t)x 

Pn 

(6/t): 

''n 

(6/t): 

X Pn 

(6/t)x 

■ 

(6/t); 


100 

100 

100 

IOC 

100 

0.02 

33.3 

0.0 2 

33 . 3 

0.02 

33.8 

0.02 

34.1 

0.03 

34.3 

0.11 

109.8 

0,12 

127.9 

0.13 

127.4 

0.15 

125.5 

0.16 

123 .9 

0.57 

220.4 

0 . 63 

205 .0 

0.67 

188.3 

0.75 

172.5 

0.84 

166.4 

2.89 

288.4 

3.18 

272 .2 

3.39 

261.2 

3.78 

242.4 

4.24 

230.2 

L4.45 

578.6 

15.93 

574.0 

16.99 

586.0 

18.91 

545.4 

21.22 

531.4 

>3.70 

667.4 

28.05 

919.6 

30.47 

737.4 

36.31 

758.6 

41.82 

764.6 


d. H 

• = 100x10 N/m^ 







c/t=0.C3 

c/t=C 

.06 

c/ 1=0 

.085 

c/t=0 

.14 

c/t=0.20 

(<5/t): 


(S/t): 


(6/t); 


(6/t): 


( 6/t)x 

■ 

100 

100 

100 ■ 

100 

100 

0.02 

33.3 

0.02 

33.3 

0.02 

33.8 

0.02 

34.1 

0.03 

34.3 

0.11 

131.2 

0.12 

129.0 

0.13 

128.5 

0.15 

126.5 

0.16 

104.8 

0.56 

225.6 

0.63 

213.4 

0.67 

195.7 

0.75 

178.3 

0.84 

172.1 

2.88 

313.2 

3.18 

297.8 

3,40 

288.2 

3.78 

269 .0 

4.24 

256.4 

L4.44 

699.4 

15.93 

699.6 

16.99 

696.4 

18.91 

672.8 

21.22 

657 .6 

24.04 

834.0 

27.79 

886.4 

30.47 

917.0 

35.91 

950.8 

41.82 

972.8 
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Table 4.2 Punch Penetration vs Punch Force Variation at 

Different Fractional Clearances (t=3xlO“^rni 6xlO“^m) 

E=12500xlo'^ N/m^, y = 0.34, H’ =25x10”^ N/m^, 8^=1. 25, 

a =QxlO^ N/m^ 
yp 


a. t = 3xl0~^ m 


c/ t=0.03 

c/t=0, 

.06 

c/t=0.085 

c/t=0.14 

c/t=0.20 


(6/t)x 

IOC 

Pn 

(6/t): 

100 


(6/t)x 

100 


(6/t)x 

100 

"n 

(5 /t)x 
ICO 

Pn 


0.02 

65,6 

0.02 

68.6 

0.02 

69.2 

0.02 

69.7 

0.03 

70.2 


0.11 

257.8 

0.12 

256.7 

0.13 

255.3 

0.15 

251.7 

0.16 

249.7 


0.59 

414.0 

0.64 

367,8 

0.68 

352.4 

0.75 

327.5 

0 .85 

313.3 


2.94 

466.4 

3.22 

441.4 

3.40 

422.2 

3.79 

390.9 

4.26 

362.9 


14,73 

665.9 

16,14 

645.6 

17.00 

628.8 

18.94 

588.1 

21.33 

565.3 


23.69 

7C2.7 

27.57 

714,2 

29,86 

707.2 

35.25 

693.4 

41.33 

695.8 


b. t 

= 6x10' 

m 

HI 


7^r>ci 


1 4 

^ - - 



c/t=0 

l6/T)x 

100 

.03 

c / t=0 .Ub 
(6/t)x Pj. 
ino 

c/ 

(6 /t}x Pj^ 
100 

c/ 

Ca/t/x 

100 

• in’ 

U/t)x 

100 



0,.02 

132.0 

0.02 

136.6 

0.02 

136,6 

0.02 

139.0 

0.03 

139.8 


0.11 

516.2 

0.12 

512,8 

0.13 

509.6 

0.15 

501.4 

0.16 

496.8 

1 

0.59 

828.6 

0.64 

772.0 

0.67 

745.2 

0.75 

688.2 

0.84 

625.3 

' 

2.96 

933.8 

3.22 

899.2 

3.38 

869.2 

3.80 

805.1 

4.24 

720.8 


14.83 

1335.8 

16.14 

1316.8 

16.92 . 

1274.0 

19.00 

1206.4 

21.21 

1120.6 


23.80 

1407.4 

26.86 

1422.2 

29.27 , 

1424.0 

34.97 

1419.4 

41.24 

1388.7 



83 


punch penetration at crack initiation is plotted against the 
variation of percentage clearances for three different values 
of H' • It is observed from the results that the value of 
percentage punch penetration at crack initiation increases with 
the increase of c/t ratio. This may be due to the fact that 
with larger clearance the material can deform easily > i.e., 
the severity of deformation decreases, and, thus, requires 
more punch penetration before crack initiation. With a very 
high value of c/t ratio more of drawing action takes place. But 
it can bo soon that for the same percentage clearance, by vary- 
ing the value of H' the percentage penetration required at 
crack initiation remains nearly constant. Similar observations 
are made from the plots of the percentage penetration at crack 
initiation vs. percentage clearance (Fig. 4.9) for different 
sheet thicknesses, i.e., by varying the sheet thickness for a 
specific value of c/t ratio the percentage punch penetration 
at crack initiation remains more or less fixed. So it can be 
concluded that the percentage punch penetration at crack ini- 
tiation is insensitive to the work hardening parameter, H , a 
the sheet thickness, t. But it strongly depends on the local 
fracture strain of the material of the sheet. The above ta 
ment becomes clear when Fig. 4.10 is examined critically, it 
can bo seen that the percentage punch penetration at crack 
initiation increases with the increase in the value of local 
fracture strain of the material. This is observed at a 


clearances. 
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4.2.3 Maximum Punch Force 

From the results of the analysis an attenpt is also made 
to study the effects of various parameters on the maximum punch 
force, (i.e., the punch force required for crack initia- 

tion). It is described here in brief with the help of graphs. 

In Fig. plotted against the variation of sheet 

thickness. As the sheet thickness is increased P„_„ increases 
linearly with it. Similar trend has been also observed from 
experimental results [25] and even in the analytical work 
reported [23]. Fig. 4.12 illustrates the effect of varying 
c/t on P„^„ for various values of H* . Here the value of P^-„ 
is found to have little variation at smaller values of the c/t 
ratio but is otherwise almost independent of c/t ratio, within 
the range of the clearances used. This observation was also 
made by a number of workers in this field, and is already men- 
tioned in Chapter 1. But maximum punch force is seen to 
increase with the increase in H’ . 

4.2.4 Optimum Clearance 

In the present work, one of the objective is to arrive at 
a suitable value of c/t called the optimum clearance (c^/t) from 
the results of the analysis. As mentioned earlier it is postu- 
lated that the direction of the crack propagation coincides with 
the direction of the maximum shear strain. When the crack ini- 
tiates the diagonal angle (9) of the deformed continuum and the 
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direction of tho maximum shear strain (0) are calculated with 
w.r. to y-axis for various c/t ratios. Tho results are tabu- 
lated in Table 4.3, and Fig. 4.13 shows one set of such results 
graphically for = 1,1. As c/t ratio is increased the angle 

0 also increases in a direct proportion but the angle 0 decrea- 
ses. Thus both tho plots intersect each other at some point 
and the value of c/t corresponding to this is taken as the 
optimum clearance ratio (c^/t), for the values of the para- 
meters used. The term optimum clearance is used here because 
when tho direction of the crack propagation coincides with 
the diagonal, the fractures meet resulting in clean edges. The 
amount of tho work required to be done for conplete separation 
of blank will also be minimum. In Fig, 4,14 and Fig. 4.15 
similar graphs are drawn from the results (Table 4.3), for the 
values of = 1,25 and = 1.59, respectively. The corres- 
ponding optimum clearance values are also found in the same 
manner mentioned above. The values of the optimum clearance, 
thus obtained, are plotted against the 8^, as shown in Fig. 
4.16. It is seen from this plot that as the ductility of the 
material (i.e., e^) is increased optimum clearance required 

between the punch and tho die reduces e.g., the optimum 
clearance required for blanking aluminium sheet (higher e^) 
is less than that required for blanking mild steel and brass 
(lower 8^). When the results obtained by the finite element 
analysis are compared with the experimental values [25] and the 
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Table 4.3 Maximum Punch Penetration, Diagonal Angle and 
Direction of Maximum Shear Strain for 
different Fractional Clearances and Ductility 

t = Sx!©”^ m, o = 8xlo”^ 

yp 

E = 12500x10*^ N/m^, y= 0.34, H* = 25xlo'^ N/m^ 

0 = Direction of maximum shear strain with y-axis 
9 - = Diagonal Angle with y-axis 


c/t 

e f ^ ^ 

.1 

Bf 

= 1.25 


= 1.59 


^100 

0° 

0° 

/A+Cs 

(-T")x 

100 

0° 

0° 

(^)x 9° 

100 

0° 

0,03 

20.42 

2.15 

13.60 

23.69 

2.25 

14.56 

29.8 2.44 

13.58 

0.06 

23.59 

4.49 

12,00 

27.57 

4.73 

12.54 

34.56 5.23 

11.80 

0.085 

26.04 

6.56 

11.41 

29.86 

6.91 

11.10 



0.14 

30.83 

11.44 

10.40 

35.25 

12.20 

10.64 

44.83 14.23 

10.30 

0.20 

35.97 

17.34 

9.40 

41.33 

18.82 

9.78 

53.05 23.07 

9.38 



theoretical curve from the previous work [23] , from the same 

figure, the trends are observed to be similar and these resuli 

match better with the experimental curve- Afterwards using i: 

plot of Fig. 4.9 the values of the optimum percentage punch 
A + c 

penetration ( (~”:^ — 2)xl00) for all the values of used are 

obtained, corresponding to the respective c^/t ratio. Fig, 4,1 

A t Cq 

shows the variation of ( ;£—*-) xlOO with the variation of 

It can be seen that the value of the optimum punch penetration 
is around 30?^, for the range of selected. When the experi 

mental [25] as well as the theoretical [23] curves are super- 
imposed on the finite element solution, it is seen that the 
results obtained by the finite element method match better wit] 
the experimental values. 

An. attempt was also made to get the crack path in the 
continuum with the help of the existing programme, which is 
described in brief here. The crack initiation criterion sta- 
ted earlier is applied after each displacement step, to chock 
whether crack is initiated at any point or not. Once the 
indication of crack initiation is obtained, the c ontribution oi 
that element for the overall stiffness of the continuum is made 
zero. This is done by assigning zero elastic modulus value 
(i.o., E = 0) to that element, thus the overall stiffness of 

I 

the continuum decreases. Then the next convenient increments 
of the displacement steps are applied and the analysis is ' 
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4.2.5 Deformation Pattern and Strain Distribution 

Nodal displacements, total stresses and total strains are 
obtained after each increment step at various points in the 
material, from the results of the analysis. From these the 
plots of the deformed grid pattern and the distribution of the 
strains in the clearance gap are drawn. All the figures 
which are described below are shown with separate scales for 
the clearance gap and the other part of the continuum. The 
clearance gap is exaggerated for clarity and better under- 
standing . 

In Fig. 4.18 to Fig. 4.22 the deformation of the various 
points, from the beginning of the punch penetration till the 
crack initiation, for c/t ratio of 0»20 is shown. It is seen 
that until the punch penetrates up to about 0.85^ into the 
sheet there is no appreciable deformation of the nodal points. 
The deformation becomes pronounced only with the further punch 
movement. It is seen that the points near the die and the punch 
corners undergo maximum deformation, while for those away from 
the corners and towards the centre of the sheet the deforma- 
tion is less severe. This brings out clearly that the crack 
should initiate at the punch and the die corners. The effect 
of change of c/t ratio on the deformation behaviour (at crack 
initiation) is brought out in Fig. 4.22 to Fig* 4.25. The 
deformation is more severe for small -clearance and becomes less 
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severe with increase in the clearance# This indicates that for 
the same value of s^, the percentage punch penetration required 
(upto crack initiation) increases as the clearance is increased. 

In Fig. 4.26 to Fig. 4.28 the variation of the maximum 
and the minimum principal strains, and the effective strains, 
in the clearance gap are shown for all the displacement steps. 
They are drawn in a part of the sheet thickness and only near 
tho die for c/t = 0.20. It can be seen that, the contours 
change abruptly till the percentage punch penetration is about 

0. 847, then there is no further significant change in the shape 
of tho contours. The patterns for the fourth, fifth and the 
sixth steps look similar in each figure. It is found that the 
highest numerical values of these strains occur near the lower- 
most loft corner, i.e., the die corner, which indicates that 
the crack should initiate from that corner only. The gradient 
of the contours are steepest at tho die corner and decreases 

as one goes away from it. Fig, 4.29 to Fig, 4.31 show 
corresponding strain contours for four different c/t ratios. 

Tho maximum value of the strains are observed near the die 
corner, for any c/t ratio. But it can be seen that strain 
distribution in tho clearance gap is affected by the c/t ratio, 
and for c/t ratio in the range of 0.14 to 0.20, the lowest 
strain occur near the lower right portion of the clearance gap, 

1. e., opposite to the die comer. 
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4,2.6 Spread of the Plastic Zone 

In order to study the spread of the plastic zone in the 
sheet material, in relation to the punch penetration, an analy- 
sis is performed by prescribing small displacement steps. After 
each step the number of sampling points which yield are found 
out. The region enclosed within the points gives the shape 
and the size of the plastic zone in the material. The plastic 
enclaves so obtained are superimposed one above the other as 
shown in Fig. 4.32. It shows that the plastic zone expands 
slowly in the beginning, and when the value of the percentage 
punch penetration is about 0.11 all the material in the 
clearance gap is almost completely plastic. This value of 
penetration is very low when compared with the value of per- 
centage punch penetration required for crack initiation (which 
is of the order of 41.335^). Afterwards the plastic deforma- 
tion of the material in the clearance gap continuously 
increases with further punch movements. This indicates that 
during the process of blanking severe plastic deformations are 
encountered in the clearance gap. 
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Fig. 4.2 Variation of punch force with punch penetration 
for t = l.47XI0’m and H' = ZSXlo’^N/m^ 



94 


Fig. 4.3 Variation of punch force with punch penetration 
for t = 1.47 X lO'^m and H' = 50 X lO^'N/m^ 


(6/t)XI00 

Fig. 4.4 Variation of punch force with punch penetration 
for t = I.47X la’m and H'= 75 X 10^ N/m^ 
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Fig. 4.6 Variation of punch force with punch penetration 
for t = 3XI0^m and H'=25XI0'^ N/m^ 






(6/t)X100 


Fig. 4. 7 Variation of punch force with punch penetration 
for t = 6Xl6^m and H'=25X10^ N/m^ 
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Fig. 4.8 Variation of punch penetration (at crack 
initiation) with percentage clearance. 





Fig. 4.9 Variation of punch penetration (at crack 
initiation ) with percentage clearance. 
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Fig. 4. 10 Variation of punch penetration (at crack 
initiation) with local fracture strain. 



max 


102 



tm 


Fig. 4.1 1 Variation of maximum punch force with 
sheet thickness. 
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Fig. 4. 12 Variation of maximum punch force with percentage 
clearance. 
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Fig. 4.13 Variation of diagonal angle and direction 
of nnaximum shear strain with percentage 
clearance for €f = 1.1 
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Fig, 4.15 Variation of diagonal angle and direction of 
maximum shear strain with percentage clearence 
for €. = 1,59 
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Fig.4.16 Variation of optimum clearance with local 
fracture strain 
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Fig. 4.17 Variation of optimum punch penetration with 
Ideal fracture strain. 
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Fig.4.20 Deformation pattern at (j)X 100 = 4.26 for c/t = 0.20 
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Fig.4.2l Deformation pattern at (-f-) X 100 = 21.33 for c/t = 0.20 



113 



Ffg.4.22 D#forrmitlofi pattern at crack Initiation for c/t*0.20 



t=0.003m,H' = 25XI0’'N/m2,ef = l.25 
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Flg.4.23 Deformation pattern at crack Initiation for c/t = O.I4 



= 0.003m, H’= 25 XI0’'N/m2,6f = 1.25 
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Fig.4.24 Deformation pattern at crack initiation for 
c/t:=0.06 
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Flg.4.25 Deformation pattern ot crock initiation 
for c/t* 0.03 
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Flg.4.23 Deformation pattern at crack initiation for c/t = 0,(4 
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CHAPTER 5 


CONCLUSIONS AND FUTURE \VORK 

5.1 INTRODUCTION 

As mentioned earlier the main objective of this study was 
to evolve a finite element model for analysing the deformation 
process during blanking. This model helped in determining the 
inter-dependance of various process parameters, magnitude of 
blanking load, optimum punch-die clearance and percentage punch 
penetration at crack initiation. Besides this, the plots of 
the deformation pattern, strain distribution and the growth 
of the plastic region helped in understanding the process in 
detail . 

5.2 CONCLUSIONS 

On the basis of the foregoing analyses and discussions 
the following important conclusions can be drawn : 

(i) The developed elasto— plastic finite clement model yields 
acceptable results upto crack initiation. 

(ii) By studying the growth of the plastic region with pro- 
gress of the punch it can bo concluded that the material 
within the clearance gap becomes plastic with a very small 
punch penetration (compared to the punch penetration 
required at crack initiation). 
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(iii) Theoretically obtained punch load agrees well with the 
results when applied to specific cases studied earlier. 

(iv) The maximum punch force required increases with the 
increase in sheet thickness and work hardening parameter 
but is independent of the clearance within the range of 
clearances used. 

(v) The percentage punch penetration at crack initiation 
increases with increase in punch-die clearance and 
ductility of the material, but the variation of sheet 
thickness and work hardening parameter has very little 
effect on it. 

(vi) The percentage punch penetration at crack initiation 

, for optimum punch-die clearance does not depend much on 
the parameters like sheet thickness, work hardening 
parameter and ductility, and it remains around 30 % of 
the sheet thickness. 

(vii) The optimum punch-die clearance depends on the local 
fracture strain of the material. 

(viii) The results yielded by the present model are closer to 
the experimental results (obtained in previous resear- 
ches) in comparison to the results from the previous 
theoretical work. 
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(ix) The mode of the metal deformation and variation of strains 
within the sheet, thickness, as seen from the distribution 
Patterns, are found to be most severe at the punch and the 
die corners, facilitating the initiation of cracks at 
these two points* 

(x) The hypotheses regarding crack initiation and crack propa- 
gation direction, proposed in the work, yield consistantly 
acceptable results. So these hypotheses may be considered 
acceptable. 

5.3 PRACTICAL UTILITY OF THE MODEL 

Besides helping in understanding the deformation process 

this model has the following practical utility : 

(i) Theoretical analysis, using this model, can give punch 
load and stress distribution in the sheet metal, wtiich 
can be used for selecting proper die and punch material. 

(ii) Prediction of proper amount of punch-die clearance (an 
important parameter of die design) is possible using 
this analysis, if the sheet thickness and its material 
properties are specified. Otherwise to obtain an optimum 
punch-die clearance a good number of actual experiments 
are required to be performed and their results required 
to be analysed in detail. 
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5.4 RECX)MMENDATI0NS FOR FUTURE WORK 

The work, presented in the thesis has been conducted 
using a number of simplifying assumptions but for the future 
work it is recommended to perform similar type of analysis, 
taking into account the actual operating conditions during 
the process, as mentioned below. A closed contour blanking 
process should be analysed with the present technique 'using 
three dimensional analysis, and the effects of the friction, 
tool shear, dullness of cutting edg^ clamping forces etc. should 
be considered. Anisotropy introduced in the sheet material 
during the process, the speed of the blanking and the strain 
rate effects should also be included in the finite element 
model so as to obtain the more accurate analysis of the 
actual operation. By using a suitable element and a very 
fine mesh, attempt should be made to get the complete crack 
path, i.e,, till the blank is separated from the sheet. 
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Appendix A 


ELASTO-PLASTIC PROGRAMME 

This computer programme is developed on DEC-1090 computer 
system. It is a general purpose programme in which the 'initial 
stress method' is used for nonlinear elasto-plastic analysis. 
This incorporates von-Mises yield criterion for two dimensional 
plane strain problems. The continuum is discretised using 
eight noded isoparametric element. 

The various sub- routines of the programme carry out most 
of the computational steps and rest of the calculations are done 
in main programme. Main programme controls various subroutines. 
Flow chart for MAIN PROGRAIvlME and STRESS ROUTINE are given in 
Fig. A.l and Fig. A.2 respectively. 

MAIN 

PROGRAMJ\4E 

Calls GDATA 

initialises total displacement, strain, 

stress and nodal load vector 

calls LOAD 

calls FORMK 

calls SOLVE 

calls DIS 

calls STRESS 

prints total displacement, total strain, total stress 

and not nodal load once solution converges 
c'hecks for crack initiation at integration points. 
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GDATA 

ROUTINE 


LOAD 

ROUTINE 


STIFT2 

ROUTINE 


reads and prints number of nodes, number of 
elements, number of nodes per element, order of 
numerical integration, number of nodes on which 
geometric boundary condition are specified, maxi- 
mum number of iterations in an increment, number 
of different materials, bandwidth of the global 
stiffness matrix 

reads and prints material properties 
reads and prints nodal point data 
reads and prints element data 
reads and prints boundary data 

initialises load vector 

reads and prints applied external 

load, if any, assembles load vector 

initialises element stiffness matrix 

calls GAUSS 

calls QUADS 

calls F0R^U2 

calls DTXY 

forms strain-displacement matrix 
forms stress-strain matrix 
forms element stiffness matrix 
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QUADS 

ROUTINE 


FORMK 

ROUTINE 


SOLVE 

ROUTINE 


DIS 

ROUTINE 


STRESS 

ROUTINE 


calculates shape function 

calculates x and y derivatives of shape function 

initialise global stiffness matrix 
calls STIFT2 

assembles global stiffness matrix 

applies prescribed geometric boundary conditions 
at nodes 

carries out the triangularisation of the banded 
symmetric matrix by Cholesky deconposition 
solves for displacement vector by back substitution 
for all iterations 

prints incremental nodal displacements (if required) 
calculates incremental strains 

initialize residual load vector, 
calculates incremental stresses, 
checks for yielding at integration points, 
calculates elasto-plastic stress-strain 
matrix if integrating point yielded, 
computes the correction stress vector 
and find residual load vector. 



6DATA 


JNITIALISE TOTAL STRESS, STRAIN. DISPLA- 
CEMENT AND NODAL LOAD VECTORS. 


LOAD 


INCR - I 


FORMK 


ITER - I 


SOLVE 


STRESS 


CONVERGENCY CRITERION 


ITER « MAX ITER ? 


ITER- ITER -f I 


PRINT TOTAL DISPLACEMENTS, 
STRAINS.STRESSES, NODAL LOADS. 



"tRACK INITIATIONS 
.CRITERIA SATISFIED 


INCR-INCR + 



PRINT CRACK INITIATED, GAUSS POINT 
NUMBER, ELEMENT NUMBER. 



FI6.A.I FLOW CHART; MAIN PROGRAMME 













